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RADIATIVE BEHAVIOR OF A GAS LAYER SEEDED WITH SOOT 

by Robert Siegel 
Lewis Research Center 


SUMMARY 

The heat transfer was examined for a gaseous layer seeded with radiation-absorbing 
carbon particles and flowing along a surface- The layer is subjected to an external 
high-temperature source of blackbody radiation. The particle absorption coefficient was 
taken to be inversely proportional to wavelength, which is approximately true for soot at 
wavelengths below about 5 micrometers. The radiative behavior was foimd to depend on 
a parameter containing particle concentration, layer thickness, and source temperature. 
Only a very small particle volume concentration, in the range of 10“^, was required to 
obtain high absorption in a 1-centimeter-thick layer for typical conditions. The results 
provide the distance along the surface for which the heat transfer to the wall remains 
within an acceptable limit and the particles remain below a temperature at which they 
will melt or vaporize. The wall protection by the layer lasts only until the particles 
vaporize or the layer becomes so hot that it reradiates substantial energy to the wall. 
Depending on the layer mass velocity the protection may be effective for a distance along 
the wall of only a few layer thicknesses. Hence, to protect greater wall lengths, it will 
be necessary to introduce the suspension through multiple slots or holes along the wall. 


INTRODUCTION 

In devices involvii^ high temperatures, significant amounts of energy can be trans- 
ferred by both thermal radiation and convection, thereby causing high heat fluxes at the 
bounding walls. In the case of convection, gaseous film or transpiration cooling may be 
used to reduce the heat flux reachii^ the wall. Such a protective film, however, is 
usually not effective for reducing radiative heat transfer as most gases are transparent 
in the temperature rai^e for which solid walls can exist. For example, hydrogen, which 
could be used for both propellant and coolant in a nuclear or laser rocket, is transparent 
at temperatures below about 5000 K. Common exceptions are carbon dioxide, water va- 
por, carbon monoxide, and sulfur dioxide, which absorb and radiate in distinct wave- 


length bands. However, in a thin layer, as would be typical for fihn cooling, even these 
gases would be quite transparent. 

A technique that has been proposed for reducii^ the external radiative heat load 
reaching a surface is to flow a gas layer alox^ the wall and seed the gas with small par- 
ticles that absorb radiation. This provides a heat shield between the heat source and the 
wall. Possible seeding materials that have been suggested are ground carbon particles 
or soot, hafnium carbide, tungsten, and aluminum oxide (ref. 1). In the seeded layer a 
small amount of the radiant energy that is absorbed by the particles raises the particle 
temperature. Another portion of the absorbed energy is transferred from the particles 
to the gas by conduction and convection, and in the present analysis it is assumed that 
this process is effective enough to maintain the gas and the particles at essentially the 
same temperature. The remaining absorbed energy is transferred away from the par- 
ticles by radiation. As the particle temperature increases, this reradiation becomes a 
significant portion of the absorbed energy. The wall will no longer be protected when 
significant absorbed energy is reradiated to the wall or when the particles become so hot 
that they begin to sublime or vaporize. 

If there are too few particles in the gas, radiation will be readily transmitted 
through the layer and the layer will not be an effective heat shield. A sufficient number 
of particles leads to good attenuation of the incident flux, but will also provide good re- 
radiation to the wall. A difficulty with effective radiative protection by a seeded layer 
is that the heat capacity per unit volume of the particulate suspension is generally small 
because of the low heat capacity of the gas. The particle volume concentration required 
to achieve significant opacity for the incident flux is usually very small, so that particles 
do not add significantly to the heat capacity of the suspension. This means that the gas 
temperature will rise quite rapidly in the flow direction unless the mass flow rates are 
high. Within a short distance the particles will either reach their vaporization temper- 
ature or there will be large reradiation to the wall. 

In the present report a simplified analysis is made to determine the particle layer 
heatii^ and wall heat flux dependency on factors such as particle concentration, suspen- 
sion mass flow rate, and spectral distribution of the incident radiation, which, for the 
conditions considered herein, is determined by the temperature of the radiative source. 
The analysis assumes the use of carbon powders or soot; in this instance scattering by 
the particles can be neglected in comparison with absorption (refs. 2 and 3). The soot 
spectral absorption coefficient is obtained from the Mie equations in conjunction with 
measurements of the optical constants of soot. 

The configuration being analyzed is shown in figure 1. A seeded gas layer of thick- 
ness L is flowing with an average velocity u along a wall at temperature T^. A ra- 
diative source at temperature T_ supplies a uniform diffuse energy flux q„ to the en- 

s s 

tire upper boundary of the layer. There is also convection from the layer to the wall. 

As a result a heat removal q^(x) is required at the wall to maintain the specified wall 
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temperature. It is assumed that the heat transfer between the particles and gas is suf- 
ficiently good so that the temperature difference between the particles and the gas can 
be neglected. As a result of the convective and radiative processes, the suspension has 
a temperature variation T^(x) in the flow direction. 

A nonflowing layer of suspended particles is considered first in order to develop the 
terms in the radiative heat balance and to obtain the magnitudes of the radiative param- 
eters that can provide effective radiative absorption of incident energy. One of the 
terms obtained is the net radiative energy absorbed by the suspension. This term is 
then combined with the convective terms to yield the energy equation for the local mean 
temperature of the layer flowing along the wall. This equation is solved by numerical 
forward integration, and the resulting temperature variation is used in another energy 
balance to obtain the heat flow into the wall. 


ANALYSIS FOR NONFLOWING SUSPENSION 

Figure 2 shows a stationary layer subjected to a radiative heat source at tempera- 
ture Tg. Obtaining the relation between the wall temperature and the heat removal at 
the wall will give the heat removal required to keep the wall temperature at an accept- 
able limitii^ value. The results will be expressed in terms of a parameter containing 
the particle concentration in the suspension. The term q^ is the energy per unit area 
of the wall that is beii^ absorbed by the suspension and then removed by some means 
other than radiation to maintain the suspension at a specified T^^. This heat removal 
would ordinarily be principally by the convective transport of the suspension. In the ab- 
sence of flow or any auxiliary means of internally cooling the suspension, the suspension 
will come to an equilibrium temperature sufficiently high so that all energy absorbed is 
radiated away. In this section, q^^^ is determined in terms of a specified T^^^. This 
relation is used later for the flowing layer to give the energy absorbed in the suspension 
in terms of its local temperature. 


Net Radiation Method 

The interaction of the radiative terms can be obtained by using the net radiation 
method (ref. 4), which relates the incomii^ and outgoing fliixes at each layer bomidary. 
Since the radiative properties of the particulate suspension depend on wavelength, a 
spectral analysis is required. The relations obtained will relate the wall heat flux and 
the energy absorbed in the suspension to the incident flux, the wall temperature, and the 
suspension temperature. A heat balance at the wall gives, for a small wavelei^th inter- 
val centered about wavelei^h X (where q^^ is energy flux per unit wavelei^th interval). 
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w *^Xi, w ~ ^Xo, w 


( 1 ) 


The outgoing radiation consists of emitted and reflected portions 

^Xo, w “ ^X, w®Xb, w ^X, w^Xi, w 

” ^X, w®Xb, w ” ^X, w^*lxi, w 

where Kirchhoff’s law has been used to substitute 1 - w ^X, w* symbols are 
defined in appendix A. ) 

The radiative quantities at the two boundaries of the layer are related by the sus- 
pension spectral transmittance The suspension spectral absorptance is equal to 
1 - and from Kirchhoff's law the suspension spectral emittance is equal to the spec- 
tral absorptance. The heat flux leaving the seeded gas layer consists of transmitted and 
emitted portions 


<1X0, ^ = *1X0, w^X + ®Xb, m^^ ■ '"x) 

Similarly, the incident energy at the wall is 

‘Ixi, w = ‘Ixi, + ®Xb, m(l - ^x) 

where q^j ^ is the incident spectral energy flux arriving at boundary I (fig. 2) from 
the soiuTce^ external to the layer. 

Another relation is a balance on the seeded layer between the incoming and outgoii^ 
spectral quantities 


^X, m ^Xo, w ^Xi, I ~ ^Xi, w *^Xo, I 

where q^^ is the net spectral energy absorbed by the suspension. 

Since’ it is desired to obtain relations between the heat flow quantities and the sus- 
pension temperature, wall temperature, and incident energy from the external source, 
equations (1) to (5) are solved to obtain q^^ q^^ and q^^^ ^ in terms of e^^^ 

e and q. • , . The algebraic manipu^tion is straightforward and yields the results 

*lx, w “ ^X, w^®Xb, m " ®Xb, w^ ^X, w^X^^^Xi, I ®Xb, m^ 
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m I ^X, w®Xb, w ” ®Xb, ^X, ~ ^X^X, w^®Xb, w *^Xi, I ~ ^®Xb, 


^x( 


, w *^Xi, 1 ~ ^®Xb, m^ 


X, w^^*^Xi, 1 ~ ®Xb, m^ 

( 7 ) 

w^^X^*lxi, 1 ~ ®Xb, m^ 

(8) 


Equations (6) to (8) all involve spectral energy quantities; that is, they are quanti- 
ties per imit of the small wavelei^th interval dX centered about a wavelength X, To 
obtain energy fluxes, an integration must be made over all wavelengths. The int^ration 
of equation (6) will now be carried out; then the integration of equations (7) and (8) will 
follow in a similar manner. Integrating equation (6) over all X gives 



dX = q, 


w 


X oO 

..-0 m 


■^Xb 




®Xb, m 


)dX 


( 9 ) 


To carry out the integrals, the spectral dependence of the wall emissivity ^(X) 
and the spectral dependence of the particulate suspension transmittance t^(X) must be 
known. Also the nature of the heat source must be known so that q^^. ^(X) can be spec- 
ified. The terms e^^^^ ^(X) and e^^^ obtained from Planck's law by insert- 

ing either the wall or the suspension temperature. Some specific assumptions will be 
made here that have the advantage of leadii^ to a closed-form analytical solution. This 
will reveal the general trends of how the heat flows are related to such quantities as wall 
emissivity and particle concentration in the gas. 

The wall is assumed to be gray; that is, ^ does not depend on wavelength, and 
the source providing the incident energy on the layer is assumed to radiate like a black- 
body at temperature Tg so that q^j ^ = e^^|j g. Equation (9) then becomes 

/ OO ^OO 

(«Xb, m - ®Xb, + ^w ^X^®Xb, s - ®Xb, m)*^ (1°) 

Equation (10) is normalized with respect to the source emission to yield 
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The integral of Planck’s law over all wavelengths yields the total blackbody radiation 



dX = oT^ (the Stefan-Boltzmann law). 


This is applied to equation (11), except 


in the final term. In this term an approximation will be made so that an anal 3 rtical ex- 
pression can be obtained for the numerator. The same approximation will be used in 
both niunerator and denominator so that good accuracy will be obtained for the ratio in 
this term. Equation (11) then becomes 


*w 


€ oT^ 

w s 




( 12 ) 


where 


F = 



The term is an average transmittance that accounts for the fact that radiation 
is passing through the particulate suspension in all angular directions. The shortest 
path leigth through the layer is L, which corresponds to a path normal to the wall. The 
path leigth approaches infinity for directions at grazing incidence to the wall. In travel- 
ing through an absorbing layer for a distance S, spectral radiation attenuates by a fac- 
tor e ^ where a^ is the spectral absorption coefficient. Hence the transmittance 
^ -a.S 

along the path is = e . When the integration is made for all possible path lengths 
through the layer of thickness L, the result 3 delds the average transmittance for pas- 
sage throigh a plane layer of diffuse radiation incident on one boundary of the layer 
(ref. 4, p. 563) 


~ 2Eg(a^L) 


(13) 


For all values of a^L, the exponential integral ftmction can be approximated quite well 
by the mean beam length approximation, giving 






( 14 ) 
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where Cj_^L is the mean beam length. For a plane layer the coefficient = 1. 8 is 
usually used (ref. 4, p. 570). 

If we consider now the last term of equation (12), Planck's law is given by 


®Xb = 


f C„/XT 
e ^ 



( 15 ) 


As shown for example in reference 4 (pp. 22 and 23), in the wavelength regions contain- 
ing practically all of the energy, equation (15) can be approximated quite well by Wien's 
formula 


-Xb 


2jtC 

,5 CjAT 
A e 


( 16 ) 


Substituting equations (14) and (16) into the last term of equation (12) yields 


/ oo 

F= - 


1 L 


- e 


>dK 


/ OO 


(17) 


-CgAT 
e ® dX 


Soot Absorption Coefficient 

Before equation (17) can be integrated, it is necessary to know how the absorption 
coefficient a^ varies with wavelength. This is difficult to specify in general because of 
the various t3q)es of particles that can be used, such as carbon, aluminum oxide, and 
tungsten, and the geometric forms of the particles, which can affect their optical be- 
havior. Scattering of radiation may also be significant for some types of particles, such 
as reflecting metallic particles. This would require modifying the present analysis, 
which only accounts for absorption by the particles. One approach for obtaining a^^ is 
to assume that the suspension is made of small absorbing spheres. (Soot particles have 
been found to be generally spherical, ref. 3. ) Then the Mie equations (ref. 4) yield the 
absorption coefficient in terms of the volume concentration of particles and the optical 
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constants as 


_ 3^ nK (18) 

^ (n^ - + 2)^ + 4n^K^ 

The optical constants are the ordinary index of refraction n and the extinction coeffi- 
cient K, which is a measure of radiative attenuation within the particle material. 

The calculations are carried out for soot particles as there is a reasonable amount 
of information available for soot carbon. For soot suspensions, scattering is very small 
compared with absorption (refs. 2 and 3). Also soot absorption behavior is insensitive 
to the temperature of the soot particles (ref. 3) and, hence, room -temperature informa- 
tion can be applied at elevated temperatures. In table I some measured optical con- 
stants are given for two types of soot carbon (ref. 2). To make specimens for these 
optical measurements, the particles were compressed on a flat plate. Table I also gives 
the a^/fy computed from equation (18), and these values are plotted in figure 3 

as a function of wavelength. 

The peak in the blackbody emission curve occurs at a wavelengtii (call it 
given by " ^3^^’ where Cg = 2898 micrometer K. Thus, at 1000 K the peak of 

the blackbody curve is at 2. 898 micrometers and at 10 000 K it shifts to 0. 2898 microm- 
eter. Hence, for a high-temperature system the radiant energy is at the shorter end 
of the wavelength range shown in figure 3. Shown on the figure are two lines with a 
slope of -1; and it is evident that, in the shorter wavelength region of interest, the lines 
approximate the data reasonably well. Thus, the approximate expression 


\ k 


(19) 


can be used. This is also well substantiated by the data in reference 5 and wiU be shown 
later to yield good values for the total emittance of a soot-seeded layer. The coeffi- 
cient k is found as the value of a^^/f^ when X = 1 micrometer; and figure 3 yields 

k « 4. 9 and k . , « 4. 0. For the soot from coal flames, k has been 

propane acetylene 

found to range from 3. 7 to 7. 5; for the soot from oil flames, a value of 6. 3 was found 
(ref. 6). From equation (19) the a^ depends on the product f^; the exact value of k 
is thus significant only with regard to specifying the required particle concentration. 

The form of the results obtained in the analysis does depend on the dependence on X 
as given by equation (19). For some types of particles the dependence in equation (19) 
is only a rough approximation, as indicated by the correlations discussed in reference 4 
(p. 715). 
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Integration to Obtain Final Relations 


Substituting equation (19) into (17) yields 


F = 


/ 



( 20 ) 


As shown in appendix B, this can be simplified to the form 

4 


F = 


ml 




4 \ T 


( 21 ) 


Tg + 11 


s/ 


+ 1 
m ' 


Substituting this expression into equation (12) yields the heat flux transferred into the 
wall: 


^w / m 


wi 


m 


€ 

w s 


S/ \^Sy 


^f^C^kL 


4 \ T, 


( 22 ) 


Tg + ll 


^f^CLkL 




Returning now to equation (7) and mailing the assmnptions of a gray wall and black- 
body external source give, after integrating over all wavelengths, 


/ OO y^OO y^O 


e-v dA. 
Ab, m 


- e 


w 


^ ^A^®Ab, w ®Ab, s “ ^®Ab, m^'^ “ (1 " ^ ^ ^®Ab, s “ ®Ab, m^^ 


Nondimensionalize by dividing through by / e-vy, „ dA to yield 

Jq ® 
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- e 


r 


’■x(®Xb,w + ®Xb.s-2«xb,m)^ 


W 


(1 - e^) 




dA. 


r -2 

^OO 

/ e-vu ^ 
^ \b, s '- 


Insert Wien’s law for the and equation (14) for the with the a^ from equa- 
tion (19) to obtain 


n /T \ ^ /T \ ^ 

— I I 

A ''^iT# '^iT/ 

aTg V^s/ \^s/ 


- € 


/ oo 


-d/X)[f^kCLlH-(C2/T^)] ^ ^-(lA)[f^kCLL+(C2/T3)] _ ^^-(lA)[f/CLL+(C2/Tj^)] | 


+ e 


w 


/ oo 


-CpATg 
e ^ ® dX 


dX 


- 


/ 

•A) 


j r -dA)[2f^kCLL+(C2/Tg)] -(lA)[2f^kCLL+(C2/Tj^)] 


/ OO 


- e 


e ® dX 


dX 


Applyii^ the reduction in appendix B then provides the final expression 
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In a similar fashion, equation (8) for the energy coming out of the particulate layer 
becomes 



Soot Total Emittance 


To further examine whether the assximptions used in the absorption coefficient are 
reasonable, a comparison can be made with the results for the emittance of soot sus- 
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pensions in reference 2. To do this, start with the as obtained from equation (22), 


lllllllllllllll 


a = e o 
w 


+■ 
w 


‘m 


^fyCLkL 


T3 + I 




T ^+1 
m 


(25) 


For a black wall = 1; and for low source and wall temperatures, this reduces to 


»w 






(26) 


Equation (26) represents a condition where the only radiation is from the suspension. 

Since the flux to the black wall can be written in terms of a suspension emittance as 
4 

V “ ^m°''^m’ evident that 


^m = ^ ■ 


(27) 


^fyC^kL 


T + 1' 
m 


where the Cy L is the mean path length (mean beam length) through the gas. Figure 6 
of reference 2 gives the gas emittance for propane soot as calculated numerically with- 
out using the foregoii^ assvunptions of Wien’s law and a dependence of the absorption 
coefficient on l/\. The results for are given as a function of T^^^ and of the par- 
ticle volume fraction multiplied by the path length (i. e. , in the present notation). 

With k taken from figure 3 as 4. 9 for propane soot and Cg as 1. 4388 centimeter K 
(the second radiation constant), equation (27) also gives in terms of the suspension 
temperature and the product of concentration and path length. Some comparisons are 
given in table n. The approximate emittances are in good agreement; hence, the ap- 
proximations made appear to be satisfactory to provide reasonable results for the heat 
flow quantities in the present situation. 

As can be shown from equation (27) the particulate suspension is somewhat wave- 
length selective in its behavior, meaning that its emittance for radiation at its equilib- 
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rium temperature is different from its absorptance for radiation from the incident 
source at a higher temperature. This arises from the fact that the absorption coeffi- 
cient a^^ varies appreciably with wavelength and the source radiation is predominatly 
in a shorter wavelength region than the layer reradiation. The parameter M is defined 
as M = fyCj^kLTg/C2 and hence the emittance in equation (27) is 

= 1 - l/[M(Tj^/Tg) + l] . If the suspension were at temperature T^, the emitted 
^ energy spectrum would then be the same as the incident spectrum. The emissivity 
would then become equal to the absorptance as both of these quantities would have the 
^ same weighting of their spectral properties over the incident or emitted spectrum. 

Thus when T^^ = T^ the absorptance of the layer becomes equal to its emittance and 
hence = 1 - 1/(M + 1) . Figxire 4 shows how the suspension emittance increases as 

its temperature rises. The e„ reaches the value of when T /T =1. Thus 

ill m m s 

initially when the suspension is introduced at low temperature, it will absorb well but 
will reemit poorly toward the wall. This selective behavior will aid in protectii^ the 
wall from the incident flux. 


Behavior When Wall and Seeded Gas Temperatures are 
Much Less Than Source Temperature 

To obtain some insight into the heat transfer behavior of the particle-seeded layer, 
a simplified case will be considered. If the source temperature is high, the material 
limitations of the wall and the particles will result in the ratios (T^^/T and (T /T 
being much less than unity. In this instance, equations (22) to (24) reduce to 


^w 


"W 


oT: 


^f^C^kL 


Tg + ll 


q 


m 


OT, 


= 1 


w 


1 - e 


w 




Tg + ll 


fvCykL 

\2 ~y ^ T + 1 ] 

/I o 


(28) 


(29) 
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1-6 


(30) 



w 


f„CjkL 

[2 ^ ^ T_ + 11 

/-I O 


For a black wall (e = 1), there is not any means for energy to be directed back toward 
the source. The black wall is nonreflectii^, and scattering by the particulate suspension 
has been neglected; therefore, ^ = 0. The dimensionless heat flows are a function of 
two parameters, and (f^CLkL/'C 2 )Tg; for convenience the latter parameter has been 
called M. Equations (28) to (30) are plotted in figures 5(a) to (c), and from energy con- 
servation the results in these three f^ures total unity at each M and e^. 


RESULTS AND DISCUSSION FOR NONFLOWING SUSPENSION 

Since the particulate suspension and wall are both at low temperature, figure 5 
represents the effect of radiation from the source only. Figure 5(a), which is a plot of 
equation (28), represents the fraction of energy transmitted through the particulate layer 
and absorbed in the wall. The absorption by the wall depends directly on the wall ab- 
sorptivity, which is equal to the wall emissivity for the case of a gray wall. The ab- 
scissa has in it the thickness of the layer, the particle volumetric concentration, and the 
coefficient k that is in the absorption coefficient. The figure shows how the energy 
penetratii^ throi^h the layer decreases as any of these quantities is increased. Increas- 
ing the source temperature also decreases the fraction of energy that penetrates through 
the particulate layer. This results from the absorption coefficient increasing as the 
wavelength of the radiation decreases, as indicated by equation (19). According to 
Wien’s displacement law the radiation from the blackbody source shifts to shorter wave- 
lengths as the source temperature increases; thus with increasing Tg the radiation 
shifts into a region of increased absorption by the particulate layer. 

The energy coming back out of the layer (eq. (30) and fig. 5(c)) is that which has 
reached the wall, been reflected, and then transmitted back through the layer. This 
radiation has undergone two passes through the layer, which accoimts for the factor of 2 
multiply ing the quantity containing the layer thickness in equation (30). Figure 5(b) 
shows how the absorption by the suspension increases with the parameter M. The ab- 
sorption also increases as decreases as this corresponds to increased reflection 
from the wall and thus an increase in the amount of energy that is undergoing two passes 
through the particulate suspension. 

In figure 5(a) for low M (e.g. , M < 0. 01) the suspension is becoming essentially^ 
transparent. Hence the curves giving the fraction of incident energy absorbed 
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each approach their value of e_ as M becomes very small. (Note that e is equal 
to the wall absorptivity because the wall is gray. ) When 0. 2, the energy absorbed 
by the wall is reduced to about e^/2; when M « 0.8, the energy absorbed is reduced to 
about €_/10. 


ANALYSIS FOR FLOWING SUSPENSION 

The behavior of a particulate suspension flowing along a wall will now be considered 
(fig. 6). Convective energy transport is assmned to dominate along the flow direction so 
that radiative transport can be neglected in that direction. The previous radiative anal- 
ysis can then be applied locally along the layer. Actually, as found later, the tempera- 
ture can chaise quite rapidly in the x-direction because of the low heat capacity of the 
layer, so that there can be significant radiation in the flow direction. Nevertheless the 
solution given here should indicate the length along the wall over which the particulate 
layer is effective in shielding the wall from the external radiation. The convective flux 

g from the layer to the wall is expressed in terms of a convective heat transfer coef- 
ficient h(x), which is assumed to be available from known correlations. If, for any 
particular numerical case, the suspension temperature goes above the temperature at 
which the particles can exist, the results will no longer be valid as the analysis assumes 
that the particles are present. The results can then be regarded as providing the loca- 
tion at which the suspension will cease to protect the wall. 

From equation (5), is the net amount of radiation that is absorbed in the suspen- 
sion, and an expression for is given by equation (23). The heat balance using the 
terms in figure 6 gives the net radiation absorbed as being equal to the convective terms 

‘ A“p. g * p]V } * V - V (31) 

Although the mass flow of the constituents remains constant, there are thermal expan- 
sion effects as the layer becomes heated. These effects are extremely complicated be- 
cause the upper boundary of the layer is a free boundary. No attempt was made to ac- 
coimt for these expansion effects. 

Substituting equation (23) into equation (31) and letting all quantities be constant ex- 
cept for Tjjj^ 3 deld the equation 
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After T (x) is obtained from this relation, the heat flow to the wall is obtained by 
using the radiative portion given by equation (25) in combination with the convection to 
yield 
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Equations (32) and (33) are placed in dimensionless form to yield 


dt 




1 + 


2t: 


m 


(M + 1)4 (Mt^ + 1)4 (Mt^ + 1)4 


- -w) 


m 


(2M + 1)4 (2Mtjj^ + 1)4 


- H(*m - V) 


(34) 




oT 




m 


-‘4- 


m 


(M + !)• 


(Mt^ + l) 


+ «(‘m - *w> 


(35) 


Equation (34) is a nonlinear first-order differential equation and can be readily solved by 
numerical forward integration by the Runge-Kutta method starting from an initial value 
tju = tjn, o' substituted into equation (35) to yield ^(x); this is 

the energy that must be removed from the wall in order to keep the wall cooled to the 
specified value of t^. 


RESULTS AND DISCUSSION FOR FLOWING SUSPENSION 
Raises of Parameters 

Before carrying out some typical examples it is necessary to examine what the typ- 
ical raises are of the parameters M, e, t^ ^ t^, and H. From figure 5(a) it is 
found that an almost transparent layer would correspond to M w 0. 02, one that is about 
half transparent to M « 0. 2, and one that is almost opaque to M w 1. 0. Typical 
values of might be 0. 1 to 1. 0 although it may be difficult to obtain the low values 
for the materials required in a high -temperature system. The t^ ^ is the ratio of the 
particulate suspension entrance temperature to the radiative source temperature. If 
the source were at 5000 K, a t^^ ^ of 0. 1 would mean the gas enters at 500 K. For a 
high source temperature, the t^’ ^ would typically be lower, say 0. 05 to 0. 01. For a 
wall temperature of 1000 K and source temperatures of 5000 K or above, the t^ might 
typically be 0. 2 to 0. 05. The convective heat transfer coefficient from the flowii^ gas- 
particle suspension would generally be low because the heat transfer fluid is primarily 
a gas. The h might t3q)ically be 5 to 25 W/(m^)(K) for pressures near atmospheric. 
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For a T = 10 000 K this yields H of about 10“'^ to 5x10"'^. The H values can be ap- 
preciably higher in high-pressure systems. Some typical solutions were evaluated for 
these ranges by numerical integration and are discussed in the next section. 


Typical Results for Flowing Suspension 


Some of the most important quantities that can be regulated, such as particle con- 
centration and layer thickness, are contained in the parameter M = (f^CLkL/C 2 )Tg. 
Figure 7 shows the effect of various M values on the suspension temperature and heat 
flow to the surface as functions of dimensionless length along the surface. Part (a) is 
for a black wall, and part (b) is for a gray wall with = 0. 5. Recall from figure 5(a) 
that for M = 0. 02 the layer is quite transparent, for M = 0. 2 the layer is about one- 
half transmitting, and for M = 1 and M = 2 the layer is quite strongly absorbing. For 
M = 0. 02 the suspension heats slowly in the X-dlrection because of its low absorption. 
One set of results for this M value is shown in figure 7(b), where a log scale is used 
for the abscissa because of the relatively large X involved. For this low M the heat 
transfer to the wall is essentially all by absorption of radiation directly transmitted 
through the quite transparent layer. The suspension radiates relatively little even when 
it is hot; hence, the heat flux to the waU rises only slightly as the gas temperature be- 
comes high in the region of large X. 

For a high M, such as 1 or 2, the incident flux is strongly absorbed by the suspen- 
sion, so the heat load directly transmitted to the wall is very small. At small X there 
is also little radiation by the suspension because it is still cool, and thus little radiant 
energy reaches the wall in this region. As a result of the strong absorption, however, 
the suspension heats up quickly in the flow direction and then radiates to the wall, which 
rapidly diminishes the wall’s protection from excessive heating. If the suspension tem- 
perature becomes too high, the protection is lost by vaporization of the particles. Thus, 
it is possible to shield the wall only to a limited X value by use of a high-M layer. As 
shown by figure 7 the suspension temperature variation with X is not significantly al- 
tered by increasing M beyond about 1. 0 since for this large an M essentially aU the 
energy is already being absorbed by the particles. Increasing M to 2 causes the par- 
ticulate layer to be almost opaque, as evidenced by the energy transmission to the wall 
being very small at small X. At larger X the energy is reaching the wall by suspen- 
sion reradiation. 

It is worthwhile examinii^ the physical significance of the dimensionless coordinate 


X 


ot;; 


L uRl - f )p c 
Lv v'^g p,g 


Vp'^p, p^ 


Lu^pTg 
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The numerator is the radiation incident on the layer from the inlet to x. The denomi- 
nator is the maximum energy that the layer can absorb; that is, it is the energy required 
to heat the layer from zero temperature to Tg. Thus, if all the energy incident on the 
layer were absorbed and retained in the layer, a layer starting at zero temperature 
would reach the source temperature at X = 1. This is shown by the long- and short- 
dashed curve in figure 7(a). Thus, X can be approximately regarded as the number of 
such ideal heating lengths required in the actual heat transfer process. 

To better understand the physical significance of the values of the dimensionless 
quantities, it is well to consider a numerical example. When M = 1, the layer is a good 
absorber . To obtain the magnitude of a typical particle concentration, let Tg = 5000 K 
and the layer thickness be L = 1 centimeter. Then from the definition of M, the par- 
ticle volume concentration is 

= .. 1X1. 44x10-2 mK ^ 3 

1. 8X4. 9X10"^ m X 5000 K 

Thus a very small particle volinne concentration is required to obtain good absorption of 
radiation. 

To compute the thermal behavior of the flowing layer, the ^ is needed. Assume 
that the suspension is supplied at room temperature and atmospheric pressure and is 
composed of carbon particles in nitrogen. Then 

^p = “ V^g^p, g + Vp^^p, p = kg/“^Kl- 04X10^ W secAg K) 

+ (3.27X10"^)(1300 kg/m^)(l. 26x10^ W sec/kg K) 

= 1. 186X10^ + 0. 054X10^ = 1. 240x10^ W sec/m^ K 

The particles contribute only about 5 percent to the value of ^ . If the carrier gas 
were hydrogen, the result would be about the same. (Hydrogen has about 1/14 the den- 
sity but 14 times the specific heat of nitrogen, so the product pc is about the 

\ S P> g 

same . ) 

From graphs such as figure 7, most of the layer heatii^ is seen to occur when X 
is near unity. If the layer velocity is taken to be introduced at a velocity of 10 meters 
per second, the length coordinate x for which the layer has become quite heated and 
thus loses its ability to protect the wall is 
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X = = (1)(10~^ m)(10 m /sec)(l. 24x10^ W sec/m^JC) ^ q 0173 ^ 

crT^ (5. 73X10“^ W/m^ k'^)(5000^ K^) 

= 1. 73 cm 

By going to higher pressures the mass flow can be increased, thereby increasing the ef- 
fective lei^th of the film protection. This will also increase the convective heat trans- 
fer to the wall. It appears that for many practical conditions it will be required to in- 
ject the layer along the wall at intervals of the order of the thickness of the seeded 
layer. If the injection is through holes, the wall would become somewhat like a per- 
forated plate. 

In figure 8 a layer of intermediate absorbing ability (M = 0. 2) is used to illustrate 
the effect of various wall emissivities. The dimensionless suspension temperature and 
wall heat flux are given as functions of X for three different e^. When the suspension 
and wall are both cool, figure 5(a) shows that for M = 0. 2 about one-half of the energy 
from the source is transmitted to the wall. For e = 1. 0 all of this energy is absorbed 

W A 

at the wall; hence for a small X in figure 8 the curve of q,„/oT for e_ = 1 starts at 
about 0. 5, and the curves for smaller are proportionately lower as a result of re- 
flection from the wall. As the particle suspension becomes heated it begins to reradiate 
the energy that was absorbed in it. With increasing X an equilibrium is reached 
wherein all the absorbed energy is being reradiated and the suspension temperature no 
longer increases. For M = 0. 2 about one-half of the incident energy is directly ab- 
sorbed in the suspension; and if the wall is black, it absorbs the remaining half of the 
energy. When the layer is sufficiently hot to reradiate the absorbed energy, one-half 
of this energy is reradiated to the wall and the other half back toward the source. Hence, 
for M = 0. 2 and = 1, about three-fourths of the incident energy is absorbed in the 
wall as the suspension temperature becomes high. 

The parameter H contains the heat transfer coefficient from the flowing particulate 

layer to the wall. The H values are small for the conditions of interest herein, where 

the heat transfer coefficient from the gaseous layer to the wall is not very high and the 

radiation originates from a high-temperature source. A typical set of cases were cal- 

-5 -2 

culated for H from 10 to 10 , and no significant effect was found for this range 

(fig. 9), The results are dominated by radiative effects and by the energy transport 

being carried along by the flow. Most of the results given in the other figures of this 

-4 -2 

report were calculated for H = 10 but should apply for any H less than about 10 . 

Figure 10 illustrates the effect of the dimensionless wall temperature on the sus- 
pension temperature and on the heat flvtx to the wall. This parameter is generally of 
minor importance for two reasons. Compared with radiative effects, a relatively small 
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amount of energy is being transferred by convection between the wall and the particle 
suspension (as indicated by the results in fig. 9). With regard to radiation from the 
wall, since this depends approximately on the fourth power of temperatture, has to 
be an appreciable fraction of Tg before the wall radiation becomes significant relative 
to that from the source. The curves show a small decrease in heat flux to the wall as 

is increased, as a result of decreased convection to the wall and increased radia- 
tion from the wall. 

The final parameter to be discussed is the entrance temperature of the particulate 
suspension. Results for T^^ c/'^s ^ ^ given in figure 11 for a mod- 

erately absorbing suspension, M = 0. 2, and a strongly absorbing suspension, M = 1. 

As would be expected the higher entrance temperature increases the wall heating and 
the suspension temperature at each location along the surface. 


CONCLUSIONS 

The interaction between a seeded gas layer flowing along a surface and radiation 
from a blackbody source has been analyzed. The analysis yields the energy being ab- 
sorbed by the surface and the temperature rise of the flowing seeded layer. The radia- 
tive behavior of the layer was found to be governed by the parameter 
M = (f^CLkL/C 2 )Tg and the ratio of suspension temperature to source temperature 
Tm/Ts* The parameter M contains the volume fraction of particles f^, the layer 
thickness L, and the source temperature T^, as well as the ratio of mean beam length 
to layer thickness Cj^, the constant k in the absorption coefficient, and the constant 
C 2 in Planck's energy distribution. When M » 2 the layer absorbs practically all the 
radiation incident on it, and when M 0. 2 about one -half of the incident radiation is 
absorbed. The spectral absorption coefficient of the particulate suspension is approxi- 
mately proportional to the reciprocal of the wavelength of the incident radiation. This is 
what causes the M to increase with increasing Tg as the blackbody spectrum of the 
incident radiation is thereby shifted to smaller wavelengths. 

The heat transferred into the wall can be reduced by increasing the value of M so 
that direct transmission of radiation through the layer is made small. This is a work- 
able approach only if the layer temperature stays low since with an increased value of 
M the layer will tend to emit well as its temperature increases. The heat flux to the 
wall can also be reduced by making the wall a good reflector for incident radiation. Al- 
though it was not considered here, a spectrally selective wall (e. g. , a wall with a white 
ceramic coating) would be beneficial. This wall would tend to reflect away the incident 
energy from the high-temperature source and would emit well the loiter wavelei^h 
radiation characteristic of the wall temperature. The heat transferred into the wall 
would thereby be reduced. 
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To obtain good absorption in the suspension, the M value must be at least near 
rniity. The value of M depends on the product of particle volume concentration, layer 
thickness, and source temperature. For a typical layer thickness of 1 centimeter the 
required particle volume concentrations f^ were found to be quite small, being in the 
range 10"® to 5x10”® for source temperatures of 5000 to 20 000 K. 

The results provide the distance along the surface for which the wall heat transfer 
remains within an acceptable limit and the particles remain below a temperature at 
which they will melt or vaporize. If the layer becomes hot enoi^h to vaporize the par- 
ticles or radiate appreciably to the wall, the layer is no longer effective in providing 
surface protection. The effective distances are usually short, beii^ of the order of a 
few layer thicknesses as indicated by the numerical example in the discussion of the 
flowing layer. This is because the heat capacity per unit volume is small for the com- 
bination of particles and gas. This follows from the fact that the particle volume con- 
centration required to obtain good absorption is quite small and hence the particles do 
not appreciably increase the volume heat capacity of the flowii^ layer. The heat capac- 
ity remains low as is typical for a gas, and the layer will heat rapidly as radiative en- 
ergy is absorbed. The layer heat capacity per unit volume can be raised by increasii^ 
the suspension pressure, thereby increasii^ the length over which the layer is effec- 
tive. 

In a practical design it would probably be necessary to introduce the suspension at 
frequent intervals aloi^ the surface, the spacing between introduction locations depend- 
ing on the mass flow of the layer. The wall might be made in the form of a perforated 
plate with the coolant being introduced through many holes; this would be somewhat like 
a transpiration- cooled wall but with a seeded gas to absorb and carry the incident radia- 
tion away from the wail. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 14, 1976, 

506-24. 
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appendix a 


SYMBOLS 


Cl 

Cl, C2 

C3 


F 

fy 

H 

h 

k 

L 

M 

n 


q 

T 


t ,t 
w’ m 


X 

X 


a 



K 

X 


p\ 


absorption coefficient 

ratio of mean beam length to layer thickness 
constants in Planck’s spectral energy distribution 
constant in Wien’s displacement law 
specific heat 

third exponential integral function 
radiative emissive power 
quantity defined in eq. (17) 

volvime of particles per unit volume of particulate suspension 

dimensionless quantity defined as h/oT^ 

convective heat transfer coefficient 

coefficient in absorption coefficient 

thickness of seeded layer 

parameter defined as (fyCLkL/C2)Tg 

index of refraction 

energy flux; energy per unit time and area 

absolute temperature 

dimensionless quantities T /T , T /T 
^ w' s’ m' s 

average velocity in x-direction 

dimensionless quantity defined as (x/L)^oTg/u pc^^ 

length coordinate along plate from origin of seeded layer 

absorptivity 

suspension emittance 

wall emissivity 

extinction coefficient 

wavelength 

spectral reflectivity 
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p density of gas 

Pp density of particles 

quantity equal to (1 - f>gCp^ ^ + f^pCp^ p 
a Stefan- Boltzmann constant 

T transmittance 

Subscripts: 
b blackbody 

c convective 

g gas 

i incoming energy 

I outer bovindary of seeded layer 

m gas -particle suspension (mixture) 

max at maximum of blackbody curve 

o outgoing energy 

p particle 

s source of radiation 

w wall 

X at one waveler^th; spectral quantity 
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APPENDIX B 


evaluation of integrals over all wavelengths 

The integrals in equation (20) are all of the form 

/»00 

1=1 — dX (Bl) 

where K is independent of X. Let t] = K/X. Then dq = -X"^K dX = -(r/^/K)dX. Sub- 
stitute into equation (Bl) to obtain 


I 


C /*00 


dq 


Applying this to equation (20) yields 


(B2) 



F = 


1 


\ ""2 




This is the result given in equation (21). 


(B3) 
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TABLE I. - OPTICAL CONSTANTS OF ACETYLENE AND 


PROPANE SOOTS (REF. 2) 


Wavelength, 


Acetylene soot 


Propane soot 

^9 

lim 

Index 

Extinc- 

Absorption 

Index 

Extinc- 

Absorption 

of re- 

tion CO- 

coefficient 

of re- 

tion CO- 

coefficient 


fraction, 

efficient. 

per particle 

fraction, 

efficient. 

per particle 


n 


volume frac- 
tion, 

.v>-l 

n 

K 

volume frac- 
tion, 

^x/*v’ 

-1 




fim 



jim 

0. 4358 

1. 56 

0.46 

9. 37 

1. 57 

0. 46 

9.29 

.4500 

1. 56 

.48 

9.45 

1. 56 

. 50 

9. 83 

. 5500 

1. 56 

.46 

7.42 

1. 57 

.53 

8. 44 

.6500 

1. 57 

.44 

5. 96 

1. 56 

. 52 

7.07 

. 8065 

1. 57 

.46 

5. 02 

1. 57 

.49 

5. 34 

2. 5 

2.31 

1.26 

1. 97 

2.04 

1. 15 

2. 34 

3.0 

2.62 

1. 62 

1. 44 

2.21 

1.23 

1.75 

4.0 

2. 74 

1.64 

.998 

2. 38 

1. 44 

1.24 

5.0 

2. 88 

1. 82 

.747 

2. 07 

1.72 

1. 30 

6.0 

3.22 

1. 84 

.505 

2.62 

1. 67 

.727 

7.0 

3.49 

2. 17 

.270 

3.05 

1.91 

.484 

8.5 

4. 22 

3.46 

.213 

3.26 

2. 10 

.357 

10.0 

4. 80 

3. 82 

. 143 

3. 48 

2.46 

.271 


TABLE n. - COMPARISON OF APPROXIMATE 
SUSPENSION EMITTANCES WITH VALUES 
FROM REFERENCE 2 


Concentration 
path length 
product, 

^v^L^' 

cm 

Suspension 

temperature, 

T 

m» 

K 

Suspension emittance, 

Equation (27) 

Reference 2 

0. OlXlO"'* 

1000 

0.0135 

0.013 

. 10 

1000 

. 125 

. 125 

1. 00 

1000 

.690 

.64 

.01 

2000 

.0268 

.030 

. 10 

2000 

. 232 

.250 

1. 00 

2000 

. 875 

.90 



Radiative source. 



Radiative source, 

‘li.l 



Figure 2. - Radiative quantities for stationary layer of particle-seeded gas. 
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Figure 3. - Ratio of absorption coefficient to particle 
volume concentration as a function of wavelength 
for soot particles. 


M 



Ratio of suspension temperature to source temperature. 

Rgure 4. - Emittance of suspension as a function of its 
temperature. 
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Dimensionless length. X = - — ^ 

L u pcp 

(b) Wall emissivity, 0.5. 

Figure 7. - Effect of parameter M = fyCLkLT 5 /C 2 on suspension temperature and heat flow to wall. Dimensionless 
heat transfer coefficient, H, 10‘^; dimensioniess wail temperature, y 0.1; dimensionless entrance suspension 
temperature, t™ 0.1. 







Dimensionless suspension temperature t^ ■ T^/Tj; or Dimensionless suspension temperature t 

heat flow to wall, qv^,/oT^ heat flow to wall, 


Wall 



Dimensionless heat 



Figure 9. - Effect of dimensionless heat transfer coefficient on suspension temper- 
ature and heat flow to wall. Dimensionless parameter. M. 0.2; wall emissivity, 

0. dimensionless wall temperature, t,^, 0.1; dimensionless entrance sus- 
pension temperature, t^ Q. 0.1. 



Dimensionless suspension temperature. tf„ “ or Dimensionless suspension temperature. 

heat flow to wall. q^/oTj ^ heat flow to wall. q^^/oT; 


Dimensionless 



Figure 10. - Effect of dimensionless wall temperature on suspension temperature and heat flow to mII. Dimension- 
less parameter. M. 0.2; wallemissivity. £^, 0.5. dimensionless heat transfer coefficient. H. 10"^ dimensionless 
entrance suspension temperature, t^, o> 
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